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ABSTRACT 

HD 8673 hosts a massive exoplanet in a highly eccentric orbit (e=0.723). Based on two epochs of 
speckle interferometry a previous publication identified a candidate stellar companion. We observed 
HD 8673 multiple times with the 10 nr Keck II telescope, the 5 m Hale telescope, the 3.63 m AEOS 
telescope and the 1.5nr Palomar telescope in a variety of filters with the aim of confirming and 
characterizing the stellar companion. We did not detect the candidate companion, which we now 
conclude was a false detection, but we did detect a fainter companion. We collected astrometry 
and photometry of the companion on six epochs in a variety of filters. The measured differential 
photometry enabled us to determine that the companion is an early M dwarf with a mass estimate 
of 0.33-0.45 Mg. The companion has a projected separation of 10 AU, which is one of the smallest 
projected separations of an exoplanet host binary system. Based on the limited astrometry collected, 
we are able to constrain the orbit of the stellar companion to a semi-major axis of 35-60 AU, an 
eccentricity < 0.5 and an inclination of 75-85°. The stellar companion has likely strongly influenced 
the orbit of the exoplanet and quite possibly explains its high eccentricity. 

Subject headings: binaries: visual - instrumentation: adaptive optics - stars: individual (HD 8673) 
stars: late-type 


1. INTRODUCTION 


HD 8673 (HR 410 = HIP 6702 = WDS 01262+3435) 
is identified as a member of the larger Hyades Moving 
Group although it does not fit the Hyades L i-temperature 


Budge 

1988 

). There is 

of the star. 

Saffe et al. 


(2005) derived an age of 2.8_ 


+0.5 

0.7 
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Gyr from isochrones and an 

age of 8.7 Gyr from Fe/H 

abundances. 

Holmberg et al. 

(2009 

) produced a similar 

age oi 2.5 Gyr using the photometric technique. Using 

isochrones, Valenti & Fischer 

20051, computed an age of 


4.29 Gyr. This variation in age estimates is normal due 
to a number of factors (Saffe et al. 2005) and HD 8673 
is probably somewhere m age between the Hyades and 
the Sun. It can be broadly described as a nearby solar 
type star. It was spectro scopically classified as an F7V 
dBo esgaard &; F ricl 199 01) and has a distance o f 36.1+0.5 
pc (|van Lecuwcn 2007|). |Tsantaki et al.| (|2014[) computed 
that it ad a T e // of 6472+64, a mass of 1.56+0.10 Mg 
and a radius of 1.23+0.07 Rg. 

A substellar companion orbiting HD 867 3 was detected 
via radia l velocity (RV) measurements by Hartmann et 


al. (2010). HD 8673b has a minimum mass of 14.2+1.6 
Mj U p , suggested that it is probably a low-mass brown 
dwarf rather than a planet. The companion is in a 
high eccentricity (e=0.723+0.016) orbit with a period 
of 1634+17 days. The orbit has a semi-major axis of 
3.02+0.15 AU. The detection was part of an effort to de¬ 
tect sub-stellar companions to F-type stars in order to 
increase the statistics of exoplanets around stars more 
massive than the Sun. 

In an earlier paper, we desc ribed a ca ndidate stellar 
companion to HD 8673 (Mason et al. 2011). In order to 
confirm the candidate companion and determine the im¬ 
pact of the stellar companion on the orbital dynamics of 
the planetary system, we conducted a series of adaptive 
optics (AO) observations using four different telescopes. 
Those observations now lead us to conclude that our pre¬ 
vious detection was spurious. However, we detected an¬ 
other companion that appears to be a bound low-mass 
stellar companion in a century long orbit. Details of the 
observations are in Section [2} and the analysis of the 
companion is detailed in Section [3j Finally we discuss 
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Table 1 

Astrometry 


UT 

G (°) 

P(") 

Telescope 

2004.7871 

302.3T1.0 

0.31 i0.02 

AEOS 

2011.5366 

332.3T1.0 

0.310i0.005 

Keck II 

2011.6402 

329.5T1.0 

0.308i0.1 

Hale 

2011.8671 

333.4±0.57 

0.308i0.003 

Keck II 

2012.4837 

335.2±0.63 

0.308i0.003 

Keck II 

2013.7426 

339.3±1.68 

0.32 i0.02 

Hale 


the system in Section |4j 


2. OBSERVATIONS 

We observed HD 8673 with four telescopes, each 
equipped with AO. Two of the AO systems operate in 
the visible, while two operate in the near-IR. 


2.1. KECK II OBSERVATIONS 


We used the Kec k II telescope, its AO system (Wiz- 
inowich et al. 2000), and the NIRC2 instrument to 6b- 
serve the system at three epochs in 2011 and 2012. Each 
time, we collected data in multiple filters with multiple 
coadds. Since the existence of the companion was still 
questionable, the multiple filters allowed us to determine 
if the object was an astronomical object or a quasi-static 
speckle whose position is wavelength dependent. 

After reducing each image, we measured the astrome¬ 
try and photo metry usi ng the fitstars algor ithm on each 
image (ten Brummelaar et al. 1996 2000). Photomet¬ 
ric error bars were set equal to the standard deviation 
of the measurements from all the images. We computed 
the weighted mean of the astrometry from the images 
taken in different filters on the same night, with the 
weight being the number of coadds for each image. The 
error bar was set to the standard deviation of the re¬ 
sults. That covers much of the random errors, but does 
not cover the systematic errors. We expect those to be 
smaller than those of the Advanced Electro-Optical Sys¬ 
tem (AEOS) (Section p3| ) or Hale data (Section 12 2| be¬ 
cause of the smaller point spread function (PSF) arising 
from the larger aperture size, but we do not have a way 
of quantifying those errors. The astrometry is given in 
Table [I] it lists the Besselian date of the observations, 
the position angle and separation of the companion rel¬ 
ative to the primary and the telescope used to collect 
the measurements. The photometric measurements are 
listed in Table[2] The table lists the measured magnitude 
difference in each filter, the central wavelength of those 
filters, and the telescope used to take the observations. 


2.2. HALE OBSERVATIONS 

We observed HD 8673 on 2011 August 22 UT and 2013 
September 28 UT with the Palomar Observatory Hale 5 


et al. 

2013)) 

et al. 

Arr/ id 

2001 


201 3) and the PHARO near-IR camera (Hayward 


HD 8673 was centered in the detector’s 


star in Br"f filter and after reducing the data, the frames 
were coadded. In 2013, we collected 50 frames in the K s 
filter. After the individual frames were calibrated, we 
created five images by coadding 10 frames into each im¬ 
age, allowing us to analyze the precision of the measure¬ 
ments. The fitstars algorithm was used to measure the 



Figure 1. Br'y image of the HD 8673 binary system taken with 
NIRC2 at the Keck II telescope on 12 June 2012. North is up in 
the image and East is to the left. This is a sub-image of the full 
image and is approximately l'/2 across. 


astrometry and photometry of the objects. Photometric 
error bars we re assigne d using the technique described 
in Roberts et al. (20051. The resulting; astrometry and 
photometry are listed in Tables [T] and [2] Figure [2] shows 
an image from the 2013 data collection. 



Figure 2. Ks image of the HD 8673 binary system taken with 
PHARO an “d the PALM 3000 AO system at the Hale telescope 
on 2013 September 28 UT. The companion is the object to the 
upper right. There is a ghost image caused by the neutral density 
filter to the lower left. There is a large speckle to the right of the 
primary and a mirror reflection of that speckle to the left of the 
primary. The image has the same orientation as Fig. ^ This is 
a sub-image of the full 25” field of view and is approximately 27 5 
across. 


2.3. AEOS OBSERVATIONS 
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Following the Keck II and Hale resolutions, archival 
observatio ns o btai ned with the AEOS 3.67 m telescope 
(Vigil et al. 1996) were re-examined for a comparable 
companion detection. On 2004 October 14 UT the star 
was observed i n an e ff ort to detect the bright companion 
detected by Mason et al. (2011). The observations used 
the telescope's AO system an d the Visible Imager cam¬ 
era (Roberts & Neyman 20021. We collected 1000 frames 
of the star with the Bessel 1 filter. After debiasing, dark 
subtraction, and flat fielding, the resulting frames were 
coadded with each frame weighted by its peak pixel to 
emphasize the best quality data. The result is shown in 
Figure [3j The companion is harder to see than in the 
Keck II or Hale images, due to a lower image quality be¬ 
cause of the shorter wavelength and because the system 
has a larger dynamic range in the / filter. 

Again, we measured the astrometry and photometry 
with the iterative blind deconvolution algorithm, fitstars. 
Photometri c erro r bars were assig ned using the technique 
described in Roberts et al. (2005). Tables jlj and [5] list the 
resulting astrometry and photometry measurements. 


•° 


Figure 3. I-band image of HD 8673 collected with the AEOS 
telescope and AO system on 14 October 2004 UT. The companion 
is circled. The halo of the PSF has many more speckles than the 
near-IR images, but that is to be expected due to the shorter wave¬ 
length of the AEOS images. The image has the same orientation 
as Figs. [l1and[2] This is a sub-image of the full 10 // image and is 
roughly lv8 wide. 


2.4. PALOMAR 1.5M OBSERVATIONS 


In 2012 and 2013, we observed HD 8673 four tim es with 


the Robo-AO system ( 

Baranec et al. 

2013 


2014 

) on the 

Palomar 1.5m telescope. On 2012 July 17 

UT and 2012 


July 18 UT, we used the SDSS r\ i\ z’ filters, on 2012 
October 6 UT we used the SDSS i’ & g’ filters and on 
2013 January 20 UT we used the SDSS r\ i’, z’ filters. In 
2012, the observations were collected as part of the au¬ 
tomated Robo-AO target queue, where the telescope is 
pointed at the object robotically with a 44" wide square 
field of view. For the 2013 January observations, the star 


was acquired manually which allowed the star to be po¬ 
sitioned at the exact center of camera’s field of view, and 
a 5/5 sub window of the image to be read out at 50.5Hz, 
which is 5 times faster than the standard rate for the 
full image. This improves the image quality by reducing 
the uncompensated intra-exposure stellar image motion. 
Capturing data in this manner is done infrequently, as 
it requires human intervention and drastically slows the 
acquisition rate. 

Each data set was reduced using the standard Robo- 
AO lu cky imaging pipeline to perform th e image align¬ 
ment (Terziev et a l.|2 013 


Law et al. 2014). The pipeline 


produces tour images, with each image selecting either 
the best 1%, 10%, 50% or 100% of the images. Each im¬ 
age was examined visually for companions and no com¬ 
panions were detected. We com puted the dynamic range 
of the image via the method of Turner et al. (20081 and 
determined that the images have insufficient dynamic 
range to detect the companion. We also examined the 
wider field of view image taken in 2012 and were able 
to rule out any additional companions to a limit of 10 
magnitudes fainter than the primary. 

3. ANALYSIS 
3.1. PROPER MOTION 

Between 2011 and 2012 the measurements from Keck 
II show the companion moved at a rate of 16.47mas/yr 
compared to the measu red proper motion of HD 8673 of 
251.8 mas/yr ( van Leeuwen|200 7). The low relative mo¬ 
tion of the companion indicates that it shares common 
proper motion with the primary star and is not a back¬ 
ground object. The Keck II data were used to measure 
the relative proper motion since it is a collection of data 
taken with the same instrument. The AEOS and Hale 
data are also consistent with it being a bound object. 

There is very little difference in the short term proper 
motion (a =237.19± 0.33 mas/yr, 5 — — 84.64±0.25 
mas/yr, van Leeuwen 2007) and the long term proper 
motion (a = 237 .9±1.0 mas/yr, 5 = —87.2±1.1 mas/yr, 


Hpg et al. 
canTndica 


2000). Significant differences in t hese values 
e the presence of companions (Makarov & Ka-| 


plan|2005 ), but their absence is not surprising, especially 
in the case of a pair with a large differential magnitude 
and hence large mass ratio. 

3.2. COMP A RISI ON WITH PRIOR OBSERVATIONS 

The two obs ervati ons of the companion reported by 
Mason et al. (2011) were both optained with visible 
speckle interferometry using an intensified CCD cam¬ 
era. The observations were in the Stromgren y filter 
(550±24nm) and the Am from the 2007 measure is esti¬ 
mated as 2.3±0.5. The quality of the 2001 measure was 
insufficient for an estimation of the differential magni- 
tude. Using t he same speckle in terferometry ca mera as 
Mason et al. (2011), Hartkopf & Mason (2009) did not 
detect the companion m observations using the Mt. Wil¬ 
son 2.5 m telescope on 2007.8179, only a short tim e af¬ 
ter the 2007.6049 detection in Mason et al. (2011). The 
observations had a limiting resolution "of 0'.'054. This 
would have been suffic ient to detect the O'/109 separa¬ 
tion detected by Mason et al. (2011) measured just a 
few mont hs earlier. In additio n, no companions were de¬ 
tected by Ginski et al. (2012) using the 2.2m Calar Alto 
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telescope and the AstraLux lucky imaging camera in the 
SDSS i’ filter on 2011 January 14. 

Our observations using AO have a higher dynamic 
range than speckle interferometry and would have de¬ 
tected any companions detected with that technique, 
even if those companions had considerable motion. The 


ion reported in 

Mason et al. 

(2011 

. Our non-detection 

of the candidate companion 

from 

Mason et 

al. (2011) 

and the non-detection in 

Hartkopf 

& Mason 

(2009) and 

Ginski et al. ( 

20121 lead 

us to conclude that the Mason 


known cas e of false dou bles i n speckle interfero metry; see 


McAlister et al. (19931 and Tokovinin (2012) for exam¬ 
ples and further discussion. 

The companion reported in this pape r is a r edd er ob - 
ject than the purported object in Mason et al. ([20111). 
Visible speckle interferometry using intensified GOD 
cameras can only detec t companion s with differential 
magnitudes less than 3.5 (Mason 1996). This is too low to 
detect the companion reporte d in this paper and t his ex ¬ 
pl ains why it wa s not s een by Hartkopf & Mason (2009), 
or |Mason et al. (2011). The observations of Uinski et 


|al.(l 2012|) also did not have sufficient dynamic range to 
letect the companion. 


3.3. PHOTOMETRIC ANALYSIS 


Based on the distance ( van Leeuwen|2007|) to HD 8673 
AB a nd the /-band apparent magnitude (Monet et al. 
2003J, we computed the absolute magnitude of HD 8673 
B in /-band of 9.5±0.94. U sing the 2MASS appa rent 
magnitudes in J,7J, and K (Skrutskie et al. 2006), we 
also computed the absolute magnitudes of companion in 
the J,/Z and I\ filters. Our data was taken in Jc, He 
and Ks and we do not have a transformation between 
these filters into </,// and I\. Instead we use the val¬ 
ues from narrow band filters and realize that there is an 
unestimated error in the magnitudes. We compared the 
resulting absolute magnitudes to the values listed on t he 
Dartm outh Stellar Evolut ion Databasq^q ( Potter et al. 
2008| for a 3 Gyr old star flSaffe et al .|2005 ) with a Ee/H 
of -0.01 (Nordstrom et al. 2UU4). The J,H,Ks absolute 
magnitudes all tail in the same bin. This bin has a mass 
of 0.33-0.45 M 0 , T eff of 3520-3690 K, and log G of 4.94- 
4.85. The / absolute magnitude has a much larger error 
bar and corresponds to a mass range of 0.17-0.44 M 0 . 
There is some overlap between the infrared and the I 
data. This variance in / band result is not surprising. 
The / data is the most error prone of all the measure¬ 
ments. As seen in Figure [3] the secondary PSF is in the 
halo of the primary; it also nas the largest dynamic range. 


4. DISCUSSION 

While the astrometric data shown in Table 1 does not 
cover enough of the orbit to yield a uniqe orbital solution, 
it does allow us to place some constraints on the orbital 
semi-major axis, eccentricity and inclination. The re¬ 
sults are shown in Fig. [4j We detail our procedure here, 
followed by a short discussion on the implications. 

Based on the spectral class of the primary and the 
secondary, F7V and M2V, we obtain a mass sum of ~ 
1.8M 0 (masses of 1.4 and O.4M 0 respectively). Adopting 

1,1 http: //stellar. dartmouth.edu/~models 


this total mass, we then search through the parameter 
spaces of all 6 orbital elements to look for acceptable 
solutions. Due to the shortness of the arc, there is a wide 
range of orbital elements that lead to Xreduced — 1> except 
for orbital inclinations, which are strongly constrained 
to be between 75 and 85 deg from face-on. The binary 
orbit is nearly edge-on, and the companion is now lying 
nearly along our line of sight to the star. This fact can 
be simply inferred, without detailed orbital fitting, from 
the apparent angular motion of the orbit (~ 4 °/yr), the 
apparent separation (~ 0.3”, or 11AU at a distance of 36 
pc ), and the system’s total mass. 

We can pare down possible orbital solutions using two 
further constraints. The first is the lack of evidence for 
RV residuals from the binary companion in the Hart¬ 
mann et al (2010) measurements (over 6.3 yrs). The RV 
residuals, after subtracting the signals from the highly 
eccentric planet, is shown to be nearly flat. We adopt 
a conservative limit of A RV < 200 m/s over 6.3 yrs, or 
A RV < 285 m/s over 9 yrs (the astrometry time-span). 
The second, more important, constraint arises from the 
presence of the eccentric planet - the binary orbit ought 
not to destabilize the planet (a p = 3.02 AU, e p = 0.723). 
This constraint, as we show below, rules out a large swath 
of binary orbits. 

Holman & Wiegert (1999) numerically investigated the 
stability of a planet orbiting around a star that has a 
stellar companion. In their simulations, the planet was 
initially on a circular, coplanar orbit. For a binary mass 
ratio of 0.28 (the case here), and a typical binary eccen¬ 
tricity of 0.5, the binary has to have a semi-major axis 
as > 7.1 a p for stability. This corresponds to as > 22 
AU. However, our planet has a high eccentricity, and is 
possibly not coplanar with the binary orbit. To account 
for this, we perform our own numerical simulations to 
define the stability boundary and the results, for copla¬ 
nar orbits, are also shown in Fig. [4j These allow us to 
exclude all orbits with as < 35 AU. 

While the astrometry data do not exclude orbits that 
are very wide (e.g., o_b > 60AU), these are statisti¬ 
cally less probable because the chance of catching the 
companion at the small separation of 11 AU is smaller. 
Roughly, we can peg the most likely binary orbits as 
as £ [35,60]AU and es < 0.5. 

This brings us to the question of whether the binary 
companion has contributed to the planet’s high eccen¬ 
tricity. Unless the binary orbit is wide (as ~ 60 AU) 
and circular, simulations show the companion strongly 
influences the planet’s orbit. This is even more so if the 
two orbits are substantially inclined to each other. It is 
conceivable that HD 8673b was born with a low eccen¬ 
tricity, but is pumped to its current value by the binary 
companion. It is also conceivable that the high mass of 
HD 8673b (Msini = 14.2 Mj) is related to this perturba- 
tio - the large eccentricity swing of the planet may have 
allowed it to sweep through, and accrete from, a larger 
area of the protoplanetary disk. 

In conclusion, the discovery of the companion of HD 
8673 may have provided us with an opportunity to study 
the impact of binarity on planet formation and migra¬ 
tion. There is still much work to be done on the system. 
Additional astrometric observations over the next five to 
ten years are needed to pin down the orbit of the stellar 
companion. It would also be interesting to collect addi- 
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Figure 4. All binary orbits that satisfy the astrometry con¬ 
straints (Table 1) to within deduced — ^ (blue open circles) and 
^reduced — ' (red filled circles), obtained by scanning through all 
orbital elements. All these orbits are close to edge-on (inclina¬ 
tion between 75° and 85°). Only orbits that avoid RV detection 
are shown (see text). Orbits that lie above the magenta lines (for 
anti-aligned apses between the planet and the binary orbits) will 
destabilize the highly eccentric planet, HD 8673b. A non-coplanar 
situation will also move the magenta lines downward. Admissable 
binary solutions have semi-major axis ag > 35 AU. Orbits with 
very large ag are disfavoured: 
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Figure 5. Orbital fit for the acceptable solutions in Fig. Blwhere 
red ones satisfy x 2 < 1, and green ones 1 < x 2 A 2. Thebinary 
orbit has to be nearly edge-on (with inclination angles between 
75 and 85 degrees, in order to satisfy the astrometry observations 
(points with error-bars). 


tional RV data to see if there is a radial velocity accelera¬ 
tion caused by the stellar companion, This approach has 


been fruitful in discovering low ma ss companions (Crepp 


et al. 


2012 


Knutson et al. 20141. An RV acceleration 


would further constrain the orbit. An improved orbit 
will allow us to quantitatively retrace the dynamical his¬ 
tory of the system. 
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